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^ ■ Abstract 

^ . Based on the ladder Schwinger-Dyson equation, we investigate phase struc- 

ture of the gauged Yukawa model possessing a global SU{2)i x SU{2)r 
symmetry and an unbroken (vector-like) gauge symmetry. We show that 
even when we tune the squared mass of the scalar boson in the Lagrangian 
to be positive, there still exists the dynamical chiral symmetry breaking 
due to fermion pair condensate (VEV of the composite scalar) triggered by 
the strong Yukawa coupling larger than a certain critical value. We find 
a "nontrivial ultraviolet fixed line" and "renormalized trajectories" in the 
three-dimensional coupling space of the Yukawa coupling, the gauge cou- 
pling and the "hopping parameter" of the elementary scalar fields Presence 
of the gauge coupling is crucial to existence of the fixed line. Implications 
of the result for the lattice calculation and the top quark condensate model 
are also discussed. 
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1 Introduction 

As it turned out, the standard model of modern particle physics is extremely suc- 
cessful. However, the central mystery of this model is the two missing ingredients, the 
Higgs boson and the top quark: Why are they so heavy? The ever increasing experi- 
mental bound of the top quark mass is now getting closer to the weak scale 246 GeV 
(the present LEP constraint is 164GeV ± 27GeV, while the CDF bound is > 113GeV 
[1]). This seems to suggest a special role of the top quark in the electroweak sym- 
metry breaking and hence a strong connection with yet another missing ingredient, 
the Higgs boson. In contrast to the passive role of Yukawa couplings simply pick- 
ing up the already tuned Higgs vacuum expectation value (VEV) to give mass to the 
known fermions, such a strong Yukawa coupling of the top quark may affect the entire 
dynamical picture (phase structure) of the standard Higgs sector. 

This situation can be most naturally understood by the Top Quark Condensate 
Model (Top Mode Standard Model) [2] which was proposed by Miransky, Tanabashi 
and Yamawaki (MTY)[3] and by Nambu[4] independently. This model entirely re- 
places the standard Higgs doublet by the composite one formed by the strongly cou- 
pled short range dynamics (four-fermion interaction) which is responsible for the top 
quark condensate. The Higgs boson emerges as a tt bound state and hence is deeply 
connected with the top quark itself. Actually, based on the exphcit solution of the lad- 
der Schwinger-Dyson (SD) equation of the gauged Nambu-Jona-Lasinio (NJL) model 
(QCD plus four-fermion interaction), MTY[3] predicted the top quark mass to be about 
250GeV (for the Planck scale cutoff), which is just the order of the electroweak sym- 
metry breaking scale. This model was further formulated in an elegant fashion by 
Bardeen, Hill and Lindner [5] in the standard model language: They incorporated the 
composite Higgs loop effects, which turned out to reduce the above MTY value down 
to 220GeV, a somewhat smaller value but still on the order of the weak scale. Although 
the prediction appears to be substantially higher than the present experimental bound 
mentioned above, there still remains a possibility that (at least) the essential feature 
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of the top quark condensate idea may eventually survive. 

What is the origin of the top mode four-fermion interaction, then? This question was 
first addressed in a concrete manner by Kondo, Tanabashi and Yamawaki (KTY)[6] 
(See also Ref.[7]), who suggested a heavy spinless boson exchange model with the strong 
Yukawa couphng ("Yukawa-driven Top Mode Model (YTMM)") in the framework of 
the SD equation. The idea was further discussed by Clague and Ross [8] in a slightly 
different framework. One might be tempted to consider an alternative, the heavy spin 
1 boson exchange model. However, as was already pointed out [6], this kind of model 
does not give rise to the desired four-fermion interaction, (7*^^^ term in Ref.[3], which 
communicates the top quark condensate to the bottom quark mass. Hence it has no 
chance to give mass to the bottom quark without suffering from the axion problem 
(See the second paper of Ref.[3], and also Refs.[2, 9, 10]). This difficulty also apphes 
to the more recently proposed models[ll, 12] of heavy spin 1 boson exchange. 

In YTMM [6, 7, 8] the top quark condensate is triggered by the strong Yukawa in- 
teraction. Even when we start with the same Lagrangian as the standard model (with 
opposite sign of the squared scalar mass and hence without Higgs VEV), we have 
a quite different dynamical picture: The electroweak symmetry breaking takes place 
mainly due to the top quark condensate instead of the ad hoc tuned Higgs VEV. Thus 
the principal role of the "elementary" Higgs is not to break the electroweak symmetry 
through its VEV but to supply a strong attractive force between the top and anti-top 
through the Yukawa coupling so as to trigger the top quark condensate or the composite 
Higgs VEV.i 

It should be noted that YTMM [6, 7, 8] yields potentially strong higher dimensional 
operators in addition to the original four-fermion interaction [3, 5] relevant to the top 
quark condensate. Actually, it was first pointed out by Suzuki[14] that inclusion of 



^This picture is contrasted with the Nambu's bootstrap[13] which imphes identification of the 
"elementary" Higgs with the composite one, while we here distinguish between the two Higgses, i.e., 
one (elementary) with the GUT or Planck scale mass and the other (composite) with the weak scale 
mass whose VEV makes a dominant contribution to the WjZ masses. 
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such higher dimensional operators may reduce the top mass prediction of the top mode 
standard model. The role of higher dimensional operators was further clarified by 
Hasenfratz et al.[15] (see also Ref.[16]). Once we specify a possible underlying theory, 
we can estimate the effects of these higher dimensional operators on the top mass 
prediction. For example, the spin 1 boson exchange model does not seem to yield large 
coefficients for the higher dimensional operators [12]. If, on the other hand, there exist 
large effects of the higher dimensional operators in YTMM, then this model may have 
a chance to predict a smaller top mass to be consistent with even the present LEP 
experiments mentioned above. ^ 

Apart from the top quark condensate, the strongly coupled Yukawa model may be ap- 
plied to other models beyond the standard model. An interesting example is a "Heavy 
Scalar Technicolor Model" [17] in which the role of the ETC gauge bosons in the ordi- 
nary technicolor scenario is replaced by the heavy weak doublet scalar boson exchange 
which communicates the technifermion condensate to the ordinary fermion mass. This 
is actually described by the gauged Yukawa model, with the gauge interaction be- 
ing the technicolor instead of the QCD. Such a model may have a "walking/standing 
technicolor" [18] gauge coupling. When the Yukawa coupling becomes strong, the dy- 
namical feature of this model is expected to be somewhat similar to the "strong ETC 
model" [19] based on the gauged NJL model, but may provide a rather different picture 
to be testable in the future experiments. 

In order to draw a definite conclusion on the above problems, however, we need to 
solve nonperturbative dynamics of the strong Yukawa coupling. This is a very difficult 
task, however, and cannot be done at once. Actually, as the first step KTY[6] studied 
the phase structure of the pure Yukawa theory (without gauge coupling) with a global 
SU (2) L X SU {2)r symmetry within the framework of the ladder SD equation. (The ex- 
tension has also been made[7] to include the SU{2)l x C/(l)y-invariant Yukawa model. 



^Although the YTMM might then loose predictive power for the top mass itself, it might maintain 
predictability of the top to Higgs mass ratio which would be testable in the future experiments. 
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which corresponds to the reahstic cases of the top quark condensate and the heavy 
scalar technicolor models.) It was shown that without tadpole, we have a clear signal 
of the dynamical symmetry breaking; a vanishing scalar VEV and nonzero fermionic 
condensate (fermion dynamical mass M ^ 0) for the strong Yukawa coupling larger 
than the critical coupling. However, inclusion of the tadpole correlates the both or- 
der parameters and makes the concept of "dynamical symmetry breaking" somewhat 
ambiguous, even if we have a nonzero fermionic condensate for the strong coupling 
region (Actually, there exists a critical coupling.). KTY proposed a possible criterion 
for the "dynamical symmetry breaking" that the fermionic current contribution dom- 
inates the scalar one to the decay constant of the Nambu-Goldstone (NG) bosons 
in such a way that they are mostly the composite Higgs with a small admixture of the 
elementary one. Such a situation in fact turned out to be the case in a wide range of 
the parameter space in the SU{2)l x SU{2)r Yukawa model. 

In this paper we shall extend the previous analysis [6] so as to include a vector- 
like gauge coupling ("gauged Yukawa model" in the same sense as the gauged NJL 
model) and an analog of the hopping parameter of the elementary scalar field 
whose kinetic term is parameterized as ^((?^0)^. Based on the SD equation with the 
standing gauge coupling (running effect ignored) , we study the phase structure of the 
model in the three-dimensional parameter space ( Yukawa coupling, gauge coupling, 
) compared with the previous analysis in one-dimensional space of the Yukawa 
coupling [6]. 

Inclusion of the "hopping parameter" Z^ is of course vital to our analysis, since 
this is the very parameter that characterizes the deviation, though not all, from the 
gauged NJL modeP (Z^ = 0) and hence from the original top quark condensate model. 
Actually, we find that the fermion mass is lowered due to > when compared with 
the gauged NJL model. 



^For extensive study of the gauged NJL model with standing gauge coupling, see Ref.[21] and 
references cited therein. As to the model with running coupling, see Ref.[22]. 
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Inclusion of the gauge coupling is motivated by the realistic situation of the top 
quark condensate model where the QCD coupling makes a significant contribution to 
the top mass prediction [2] . Actually, the gauge coupling drastically changes the phase 
structure: We discover the "fixed line" and the "renormahzed trajectories" which can 
only be revealed by the presence of the (standing) gauge couphng. This we believe is a 
novel feature of the present model and is one of the major achievements of this paper. 
Roughly speaking, the renormahzed trajectories correspond to the boundary which 
separates the region of the "dynamical symmetry breaking" (fermionic dominance) 
from that of the "non-dynamical symmetry breaking" (fermionic non-dominance) in 
the sense of Ref . [6] mentioned above. 

The paper is organized as follows. In Section 2 we set up our machinery, the lad- 
der SD equation for the gauged Yukawa model with a standing gauge coupling and 
a global SU{2)l x SU{2)r symmetry. We include a tadpole contribution. In Section 
3 we present a new formula ("generalized Pagels-Stokar (PS) formula") for the decay 
constant F^^ of the NG boson in this model. consists of two parts one from the 
scalar current, and the other from the fermionic one, the latter being given by the PS 
formula[20] originally developed in the theory without scalar field. The very structure 
of our formula dictates that the fermion mass must be smaller in the gauged Yukawa 
model than in the gauged NJL model. Then in Section 4 numerical analysis of the 
phase structure of this model is given. We identify the line of constant M (fermion 
dynamical mass) and constant with the "renormalization-group (RG) flow". We 
discover a "fixed line" and "renormahzed trajectories" in the three-dimensional param- 
eter space of (Yukawa coupling, gauge coupling, Z^). Section 5 is the analytical study 
of the model which reproduces the numerical one in Section 4 and further provides 
the scaling relation near the phase transition points (critical surface). We find an an- 
alytical expression for the whole fixed line. Section 6 is devoted to the conclusion and 
discussions. Possible implications for the lattice studies and the top quark condensate 
model are discussed. It is suggested that the top quark mass prediction of YTMM may 
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be lower than the original top mode prediction [3, 5] based on the gauged NJL model. 



2 Schwinger-Dyson equation 

We start with the Lagrangian of the gauged Yukawa model (Yukawa model coupled 
to a vector-like SU{Nc) gauge interaction) with a global SU{2)l x SU{2)ji symmetry: 



" c / 

,2 



+ ^tr(a^$t5''$) _ !^tr($^$) - ^ (tr($t$))' , (2.1) 

where $ is a 2 x 2 matrix, $ = cr + zt^tt" with Pauh matrices T"(a = 1, 2, 3), '0 is a 
fermion doublet field, := G^T'^ia ^ l,...,Nl - I) are SU{Nc) gauge fields, and 
''pL/R '■= Pl/riP, with Pl/r '■= (1 =F 7^)/2 being chiral projection operators. We take 
ml > 0. 

It is well known that this system with negative causes the spontaneous chiral 
symmetry breaking already at the tree level and the fermion acquires its mass in a 
passive manner from the VEV of the scalar field. However, this picture is not ap- 
propriate for the strong coupling Yukawa region where the feedback from the fermion 
sector to the scalar potential is significant as stressed before. In this region we expect 
the fermion determinant plays a more active role in the spontaneous chiral symmetry 
breaking and we need to treat the dynamics in a nonperturbative manner. Such a 
nonperturbative dynamics may be investigated by using a truncated set ( "ladder" ) of 
the Schwinger-Dyson (SD) self-consistency equations as the first approximation. 

Since a concept of renormalized parameters becomes somewhat obscure in this kind of 
nonperturbative analysis due to the ambiguity of renormalization schemes, we study the 
physical quantities (e.g., the fermion dynamical mass M, the NG boson decay constant 
Ftt, etc.) directly as functions of bare parameters (e.g., m^, (jy, etc.) and the ultraviolet 
(UV) cutoff A. Once these physical quantities are calculated, we can determine the 
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scaling properties of those bare parameters ( "renormalization group (RG)" ) in such a 
way as to fix the physical quantities on the variation of the cutoff A. Of course the RG 
functions determined in this manner are not unique and depend on the choice of the 
"physical parameters" utilized in this procedure. It should be emphasized, however, 
that the qualitative feature of the phase diagram is expected not to change depending 
upon the choice of the parameters, as far as such a choice covers (at least) minimal 
set of possible relevant operators. In this paper, we choose the dynamical mass of the 
fermion M and the NG boson decay constant for such physical parameters. 

Now, we turn to the ladder SD equation of the gauged Yukawa model. It is convenient 
to rewrite the Lagrangian Eq.(2.1) into 

9 ^\ . r- , - _ , 1 1 



+^tr (a,$taM$) _ M^tr ($t$) _ Me (tr($t$))^ , (2.2) 

where we have rescaled the scalar field $ ^ ^ Z^Nc^ and defined parameters: 
gl := N^Z^g^, M| := Z^ml and := N^Z^X^. Again note that we take M| > 0, 
since we are interested in the symmetry breaking due to the strong Yukawa coupling. 
Here we take > in order to guarantee the stability of the vacuum. Z^ plays a role 
of the "hopping parameter" in the lattice formulation (See the discussion in Section 6). 
It should be noted here that our new parameterization of the gauged Yukawa model, 
Eq.(2.2), contains a redundant parameter compared with the original one Eq.(2.1). 
We shall return to this point later. It should also be emphasized that the Eq.(2.2) 
is appropriate for the comparison to the previous analysis of the gauged NJL model 
(Z^ = A<^ = 0). 

The SD equations for the VEV of the scalar field {a) and the fermion propagator 
S{p) are given by [6] 

M|(a) + X,{a{^' + TT^)) = -j^ii^i^) (2.3) 
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and 

+^^/ ^,D'^ik)S{p-k)T'^{p-k,p) 

/ (^)^^'^'^(^ - - (2.4) 

respectively, with D^^{k) being the gauge filed propagator, D"{k) and D'^^i^k) the scalar 
and the pseudoscalar propagators, V^^, and F'^" the vertex functions. It is hopeless to 
solve a full set of the SD equations because of the lack of knowledge of the scalar/gauge 
boson propagators and the vertex functions. Here we take the ladder approximation in 
the sense that we take bare vertices and bare boson propagators instead of full vertices 
and full propagators: 

KD^ik) = r^ = -i,„ (2.5a) 

^^^"^(^^ ^ Z^fc^-M/ ""' r- = -i^,i7V, (2.5b) 
NcD,^{k) = ^(5,. -(1-0^), r^ = -i^7''7;, (2.5c) 

with ^ being the gauge fixing parameter. Though there is no solid reason to justify this 
approximation, the ladder approximation becomes plausible, when the gauge coupling 
runs slowly ("walking") or does not run at all ("standing") [18], and the effects of the 
scalar/pseudoscalar boson propagators are small.^ Actually, the ladder approximation 
yields a successful phenomenology even for the QCD running coupling used in the SD 
equation. [24] 



^The effects of the scalar/pseudoscalar boson propagators ("rainbow" graph ) are actually sup- 
pressed by \/Nc compared with those of the gauge boson propagators (gauge boson "rainbow") and 
the VEV of the scalar field ("tadpole"). (See Eqs.(2.8a)-(2.8b).) The Nc ^ oo limit yields qualita- 
tively the same phase structure as A^c = 1 case (see Section 5). 
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We substitute 



(2.6) 



into Eq.(2.4) and introduce the UV cutoff for after Wick rotation. In the ladder 
approximation, we can perform the angular integration of Eq.(2.4), which yields (see, 
e.g., Ref.[23]) 



A{yfy + B{yr 

A{y) 



(2.7a) 
(2.7b) 



The integral kernels are given by 



^B{x,y) 
K^A{x,y) 



^ KB{x,y;^) + X{l + ^/3)KB{x,y;0), 



NcZ^ 

y \ h 



X N,Z, 



2 1 

.KA{x,y;^) + -^XKA{x,y;0) \ , 



(2.8a) 
(2.8b) 



where 



2 

KB{x,y;m'^) := - / 

TT Jo 



de 



X + y — 2y/xy cos 9 + m? 
2 



x + y + m'^ + J{x + y + ni?Y — Axy 



4 /""" 

KA{x,y;m'^) := - / d9 

TT Jo 



xy cos 9 sin^ 9 



TT Jo X + y — 2y/xy cos9 + mP' 
Axy 



X 



+ y + 'm? + ^J{x + y + rn?Y — Axy 



2 ' 



(2.9a) 



(2.9b) 



and x, h and A are defined as x :— —p^i y :— —q^, 



(2.10) 



with Cp ■= YjaT"'T°' — {N^ — l)/2Nc being the quadratic Casimir of the fermion 
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representation. 

The VEV of the scalar field (u) is determined from Eq.(2.3) under factorizability 
assumption, which reads [6] 

where the fermion condensate ('^'0) is given by 

At this stage, the VEV {a) , the chiral condensate and the fermion mass function 
B{x) are inter-related among each other and are responsible for the chiral symmetry 
breaking. The C(A^) or higher terms in Eq.(2.11) do not produce significant effect 
near the critical point ^ We thus disregard those effects in the following 

calculations. Then {a) is given by the "tadpole" contribution alone: 

Combining Eq.(2.7a) and Eq.(2.13), we finally obtain the closed SD equation for the 
fermion propagator: 

(2.14a) 

r^' Aiv) 
= l+/^#Mx.y)^^(^)jfL_. (2.14b) 



In this paper we take the Landau gauge = 0). The Landau gauge is required 
to be consistent with the bare vertex function approximation at least in the pure 
gauge theory due to the Ward-Takahashi identity, since the part of the kernel Ka 
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from the gauge interaction is identically zero in this gauge and then the SD equation 
reduces to A(— p^) = 1 (see, e.g., Ref.[23]). For simplicity of calculation we take a 
further approximation for the integral kernel from the Yukawa interaction such that 
ICa — 0. For this approximation to be consistent with the bare vertex approximation, 
the solution of the coupled SD equations must lead to the result: A{—p'^) = 1 for 
A{—p'^), namely, the deviation of A(— p^) from 1 must be small. Indeed this has been 
confirmed by the previous work [6] on the Yukawa model and also in the massive vector 
boson model[25]. Thus in what follows we put 

A{-p') = 1 , (2.15) 



which implies no wave function renormalization for the fermion. Then we only have to 
solve the single integral equation for the fermion mass function B{x): 



Ah f^' yB{y) 



4:11 



y + B{yy 



+ 



/ dy)CB{x,y) 
Jo 



yBjy) 
y + Biyf 



(2.16) 



This is our basic equation. 

Although Eq.(2.16) can be solved numerically, here we solve it by converting the 
integral equation into a more tractable differential equation plus boundary conditions. 
This can be done by adopting an approximation [26] for the kernel: 



KB{x,y;m'^) 0:^ 9{x - y) — + 0{y - x 



X + 



y + 



(2.17) 



where 9{x) is the step function. Then the SD equation Eq.(2.14a) is reduced to 



Bi^) - ^ [%. ''''' 



Ml Jo ''y + Biyf 



+ / dy 
Jo 



\x Z^x + M, 



^y Z^y + M' 



yBjy) 
y + Biyf ■ 
(2.18) 
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This technical simphfication does not change the quahtative structure of the SD equa- 
tion [23]. 

Eq.(2.18) is readily converted into a set of differential equations; 

^B(x) = -(- + ^77r^^^^]v(x), (2.19a) 

-^{xV{x)) = ^^r^ (2.19b) 

plus infrared(IR) and UV boundary conditions (BC); 

V{Q) = (IRBC) (2.20a) 

with the "condensation function" V{x) being given by 

V{x)^'-rdy ^^if^ . (2.21) 

Eqs.(2.19a)-(2.20b) are the equations that our analysis is actually based on. Note here 
that the fermion pair condensate is related to V{x) as 

m = -^AV(A^). (2.22) 



3 Nambu-Goldstone boson decay constant F-j^ 

As we have mentioned before, we need to calculate physical quantity other than the 
dynamical mass of the fermion to determine the "RG" flow of the bare parameters. 
Though the mass of the physical scalar boson ■mF^^^ would be the first candidate for 
this purpose, we need to solve the SD equation for the scalar boson propagator in order 
to calculate m?'^^'*, which is far beyond the scope of the present paper. We therefore 
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choose the decay constant of the NG boson F^. 

One might suspect that the F.^ is aheady calculated as the VEV of the elementary 
scalar field (a) through Eq.(2.13). This is not true in this model, however, because of 
the presence of the mixing of tt with the composite pseudoscalar field. On might also 
think that the Pagels-Stokar (PS) formula[20] of the (composite) pion decay constant 
in QCD is applicable to this problem. However, it cannot be used as it stands due to 
the mixing with the elementary pseudoscalar field. Actually, the wave function of the 
real NG boson contains the fermion composite part (O|T''0(^)'0(y)|NG) as well as the 
elementary pseudoscalar part (0|$(a;)|NG) in the gauged Yukawa model. 

In this section, we incorporate the mixing effect into the PS formula for and obtain 
a generahzed expression for ("generalized PS formula") in such a model. We find 
that square of the decay constant F^ of the diagonahzed NG state is divided into two 
parts: F^ = F^ + F^, with Fb and Ff being the bosonic part and the fermionic part, 
respectively. The fermion part is evaluated by a certain integral formula of the fermion 
mass function B{x). This result is actually in accord with the well-known result in the 
two (elementary) doublet model: the square of the decay constant F^ is expressed by 
the sum of the (VEV)^ of each doublet. 

The NG boson decay constant is defined by 

(0|J5''(a;)|NG''(g)) = -ig^F^e'«■^5"^ (3.1) 

with |NG^(g)) being the (diagonahzed) NG boson state with momentum q [q^ = 0) 
and isospin h. The Noether current J^^ is given by 

N7j. - t" 

jan ^ _^^^tj.(^a^t9M$ _ r'^^QI^^^) - (3.2) 

The calculation of F„ becomes straightforward, if we know exactly the wave function 
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of the NG boson: 



XU ■■= (0|«f(0)|NG'^(g)), (3.3a) 
S(p)x''p^^(p,q)S{p-q) := J d'xe-'P-^{0\T^I;(x)m\NG^{q)), (3.3b) 

where S{p) is the fermion propagator and — 0. The decay constant may be 
expressed in terms of the wave function as (see Fig.l): 

N 7 

F^q^b'^' = ^ig'^tr({r«,($t)};^5,^) 

- / (0-tr(7'^75y-S(p)x^^^(p, q)S{j> - q)). (3.4) 

Generally speaking, however, it is quite difficult problem to solve the wave function 
of composite particles (Bethe-Salpeter amplitude). Fortunately, this problem is sim- 
plified for the NG boson state, thanks to the Ward-Takahashi identities of the broken 
symmetry: 



9^(rjr(x)$(y)) 

- {y,W2/))}5('H^-|/), (3.5a) 

= y75(T^(2/)^(;^))^^'H^ - 2/) + (T^(2/)^(;^))y755('n^ - A- (3-5b) 
The spectral representation of (3.5a) and (3.5b) is written as 

E(0|^(0) \n{q))-^, {n{q) \ J^mO) 
n q^ -m^ 

= {y,(^)}, (3.6a) 

/ ciV-^^-^E(0|T^(y)V'(0)|n(g))^^(n(g)|Jr(0)|0) 

= S{p)I^S-\p)^^,-^^,S-\p-q))s{p-q). (3.6b) 
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Among \n"'{q)) the only state which survives in — hmit is the NG boson state 
|NG"(g)) (m^G = 0). We therefore obtain the wave functions Xp<s>i Xp^^iPiQ = 0) 
which are completely determined at zero-momentum g = by the VEV of scalar field 
and the fermion propagator, respectively: 



Xp* = -^{y,(<f)}, (3.7a) 
XP^^(P,? = 0) = -±.{llj,,S-'ip)}. (3.7b) 



Following the paper of Pagels and Stokar[20, 9], we approximate the fermionic wave 
function 

Xp^4p,q)^Xp^^(p,q-0). (3.8) 

Now, it is straightforward to calculate the decay constant F^^ by plugging the wave 
function Eqs.(3.7a)-(3.7b) into Eq.(3.4). We find 

= F^ + Fj, (3.9a) 

= N,Z^{ay, (3.9b) 

2 Nc B-(.r) - j^D\x) , ^ 

where F^ comes from the elementary scalar wave function Eq.(3.7a) and F| from the 
fermion composite wave function Eq.(3.7b). 

Remarkably enough, this generalized PS formula definitely predicts that the fermion 
mass gets lower in the gauged Yukawa model [Z^ > 0) than in the gauged NJL 
model {Z^ = 0) for the same F^. Actually it dictates Ff ^ F^ - F^ < F^ for 
F^ — NcZ^{a)'^ > 0. Now, the fermion mass is determined by the formula 
for F|, Eq.(3.9c), as an increasing function of F|, which was the essence of the 
top quark mass determination in the top quark condensate model [3, 2]. Thus we 
have (gauged Yukawa) < (gauged NJL), since F|( gauged Yukawa) < F^ — 
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F|(gauged NJL). 

4 Numerical study 

Having set up the machineries Eqs.(2.19a)-(2.20b) and Eqs.(3.9a)-(3.9c), we are now- 
ready to perform numerical analysis of the phase and "RG" structure in the gauged 
Yukawa model. As we mentioned before, we take M| > so that the chiral symmetry 
breaking can only be caused by the fermion dynamics. Because of the redundancy of 
the parameters in the Lagrangian Eq.(2.2), we may fix one parameter — A for 
positive M| without loss of generahty. 

The goal of the numerical analysis here is to determine the bare parameters {h, A, Z^) 
as functions of the ultraviolet cutoff A and the "physical quantities" (M, F^, Xr); 

/i(A,M,F,,A«), X{A,M,F^,Xr), Z^(A,M,F,,Afl), (4.1) 

with M := -B(O) being the dynamical mass of the fermion, and A^j the renormal- 
ized gauge coupling. Dimensional analysis of Eq.(4.1) shows that h{A, M, F^^, Xr) — 
h{A/M, 1, F^/M, Xr) and Z^{A, M, F^, Xr) = Z^{A/M, 1, F^/M, Xr). Thus M can be 
set unity, M = 1, without loss of generality. Since we restrict ourselves to the standing 
gauge coupling, ^ we do not consider the RG flow to the gauge coupling A direction: 
Xr — X. We also disregard the RG flow to the A^ direction, since feedback of running 
of A^ is expected to be small as we already argued in Section 2. We thus obtain the 
"RG" flows in {h, A, Z^) space (more precisely, {h, Z^) plane sliced for each A) as we 
vary A while keeping M and F^^ flxed in Eq.(4.1). 
Numerical calculations are performed as follows in our program. For given gauge 



^It is rather straightforward to incorporate running effect of the gauge coupUng as was done in the 
gauged NJL model [22]. 
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coupling A, we solve the differential equations Eqs.(2.19a)-(2.19b); 

^B{x) = + i . ]v{x), 

^V(x) - _m + ^!^ (4 2) 

dx^^""' ~ X ^x + B^{x)' ^^-^^ 

together with the differential form of Fj — Fy(A^) Eq.(3.9c); 

d p2(._ N^ xB{x)[B{x)-^^B'{x)] 

dx ' 47r2 {x + B^{x)y ^ ' 

plus the respective IRBC's; 

S(0) = 1, nO) = 0, F|(0) = 0, (4.4) 

with (/i, Z^) being trial parameters. Then we obtain the next trial parameters (/i', Z'^ 
from Eqs.(2.20b),(3.9a)-(3.9c): 

We repeat the above steps after substituting the trial parameters {h! , Z'^) into {h, Z^), 
and so on. After several iterations, we obtain the function Eq.(4.1) with sufficient 
accuracy. 

Let us now turn to the result of our numerical calculations. We here present numerical 
result only for the U{1) case, i.e., Cp = I and A^^ = 1 in Eqs.(2.10), (3.9b)-(3.9c). 
Actually, the Nc dependence is not significant for the essential feature of the phase 
diagram as we will demonstrate by the analytical study in the next section. When 
we discuss the > 1 effect, we also have to include the running effects of the gauge 
coupling. Numerical study including both effects will be given in a separate paper. 
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We first discuss A = case (pure Yukawa model with vanishing gauge couphng). 
In Fig. 2 we show the critical line in the space (h,Zrp), which separates the chiral 
symmetric phase (M = 0) ai h < he and the spontaneously broken phase (M ^ 0) at 
h > he, with he — hc{Z(j)) being a critical coupling for each Z^.^ This is confirmed by 
the analytical study in the next section. Above the critical coupling h > he, the chiral 
condensate (ipip) exhibits a characteristic behavior shown in Fig. 3, which was first 
obtained in Ref . [6] . Note that the critical Yukawa coupling he monotonically increases 
in Zff, and the dependence of he on is very small. This also agrees with the analytical 
study. 

Next we consider another limit, = 0, which corresponds to the gauged NJL model. 
Actually, hues of the equi-correlation length — A/M in the Z^ — plane are depicted 
in Fig. 4 for A^c = 1, which converge in the ^/ ^ oo limit toward the critical line (bold 
line in Fig. 4). This critical line is nothing but the well-known critical line[27] of the 
gauged NJL model written in the space {h. A): 



for A < Ac = 1/4. (This analytical expression will be derived in the next section.) 

Then in Fig. 5 the surface of the equi-^/ is depicted in the three-dimensional bare 
parameter space {h. A, Z^) for the Yukawa coupling h, the gauge coupling A and the 
"hopping parameter" of the scalar field Z^. The ^/ ^ oo limit is the critical surface 
h = he{X, Z^) separating the chiral symmetric phase {h < he) and the spontaneous 
symmetry breaking phase {h > he), which is shown by the surface with mark "x" in 
Fig. 5. Below the critical surface {h < he), there is no nontrivial solution for the SD 
equation(M = 0). 

Now to the RG flow, which is obtained by flxing the NG boson decay constant F^^ as 



^In the previous work [6] the case of = 1 was investigated by varying the scalar mass M^. It 
was shown (see Fig. 3 of Ref. [6]) that there exists a critical Yukawa coupling he and the dependence 
of he on is very small. 




(4.6) 
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well as M. Given the value of F^^, we obtain the F^r-constant surface (two-dimensional 
manifold) in the space {h, A, Z^), which is then sliced by the A-fixed plane in accord with 
our setting of the standing gauge coupling. RG flows are obtained as the intersection of 
the F^-constant surface with the A-fixed planes, see Fig. 6. Imagine the situation that 
each RG flow passes across the equi-^/ surfaces of Fig. 5. Then each flow curve may 
be considered to be parameterized with respect to the value of We may interpret 
that the crossed points of a RG flow with the equi-^/ surfaces correspond to the result 
of successive steps of the RG transformation.''' 

All the flows in each figure of Fig. 6 are for the same F-^ and for different values of 
A. Most of the flows either run to ^ oo or ^ for := A/M ^ 1, but there is 
an indication of a single flow ( "renormalized trajectory") terminating at a finite value 
of and A ("(nontrivial) fixed point") on the critical surface. We can draw a similar 
figure for a different F^-, which indicates a different fixed point and the renormalized 
trajectory (compare (a) and (b) in Fig. 6). A set of such fixed points form a line ("fixed 
fine") which is depicted by the dotted dashed fine in Fig. 6. (The precise location of 
the fixed line is determined by the analytical study in the next section where we find 
that the fixed line extends from {h, A, Z^) = (-1/20,1/4,0) to (l/4,0,oo) for = 1.) 

We may take another look at this phase diagram by Fig. 2 (A = 0) and Fig. 7 
(A > 0), where the fiows are depicted for different values of F^ in the fixed- A plane 
{h, Z^). In Fig. 7 (A > 0) there is again an indication that a single flow (renormalized 
trajectory) terminates at a finite (fixed point) on the critical surface, while all others 
run to either Z^ — > oo or Z^ — > for ^ 1. This is sharply contrasted with the pure 
Yukawa model (A = 0) whose RG fiows are shown in Fig. 2 (sohd fines). There is no 
indication of the existence of the "(nontrivial) fixed point" . This is confirmed by the 
analytical study in the next section. This also agrees quite well with the result of the 

""Our analysis based on the SD equation can determine the flow only in the spontaneous chiral 
symmetry breaking phase where the nontrivial solution (M ^ 0, i^^ ^ 0) exists. Thus the renormalized 
trajectory cannot be drawn as extended to the symmetric phase, which of course is an "artifact" of 
this framework. We may improve this situation by calculating effective potential (or mF^^^) as we did 
in the gauged NJL model [21]. 
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lattice Yukawa model (see Section 6). Thus we found that a fixed fine exists in the 
gauged Yukawa model solely due to the presence of the gauge coupling. This is our 
main result. 

The RG fiow shows that there is a "critical value" for the decay constant which 
separates the three-dimensional space (above the critical surface h > he) into two 
regions: 

(I) < Ff {Ff := Ff{A = oo)); RG flows extend to = for ^/ > 1, 

(II) F^> F^ ; RG flows run away to ^ oo . 

The surface = Ff° is the boundary of the region (I) and (II) which consists of the 
renormalized trajectories. The intersection of the critical surface and the renormalized 
trajectories forms the fixed line. 

Finally, we consider the ratio F^/Fj for the NG boson decay constant F^, where 
F^ and F'j are the bosonic part and the fermionic part of F^, respectively, as defined 
by Eqs.(3.9b)-(3.9c). In Fig. 8 the ratio F^/Fj is shown along the RG flow (sohd 
line) in the fixed-A plane and the broken line passing through the origin denotes the 
fixed-cutoff A fine where the line with the steeper slope corresponds to the smaller 
cutoff. 

In the presence of gauge coupling, Ff is finite even in the infinite cutoff limit (Fj?° < 
oo), because the solution B{—p'^) is damping in the asymptotic UV region (see Eqs.(5.22b),(5.23b)). 
The RG flow in the region (II) runs away to ^ oo as A — > oo. However, even in 
the — > oo limit, the ratio 
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/?2 p 

— = — - 1 (4 7) 

Ff Ff ^'-'^ 

takes a finite non-zero value depending on the initially specified value of F^ for the 
RG flow. As Ftt approaches to Ff from the above, the ratio decreases and flnally 
becomes zero at A = oo just for F^ = Fj°. In the region (I) including the F^ = Ff(A'^) 
surface, therefore, we can always make the contribution from the bosonic part F^ to F^ 
arbitrarily small by choosing large cutoff or small Z^. This implies that the fermionic 
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part Fj is dominant in in the small region for large A, as was already noted [6] 
in the pure Yukawa model which has actually only the region (I) [6] . 

Comparing two points on the fixed- A line, we find a larger F^ by increasing Z^. This 
implies that the contribution of the fermion dynamical mass M to F^^ gets smaller 
by increasing Z<^, since we fixed the fermion mass M = S(0) = 1 in our numerical 
calculation. Here we notice again that the gauged Yukawa model reduces to the gauged 
NJL model in the limit — 0. Therefore this result can be rephrased as follows; by 
including the Yukawa interaction due to non-zero Z^, the fermion dynamical mass M 
can be reduced for a fixed value of F^r. This is actually a direct consequence of the 
salient feature of our generalized PS formula Eqs.(3.9a)-(3.9c), as we emphasized in 
the end of Section 3. This mechanism may be applied to a scenario of the top quark 
condensate ("Yukawa-driven Top Mode Model") [6, 7, 8] to lower the mass of the top 
quark. 



5 Analytical study 

In the previous section, wc have solved the nonlinear SD equation Eq.(2.18) numeri- 
cally. We solve here the SD equation Eq.(2.18) in a analytical manner by making use of 
the hnearization (bifurcation) [28, 23], which is valid near the critical point (M <S A). 
We thus obtain here an analytical expressions of the critical surface and the fixed line. 
These results actually confirm the numerical calculations in the previous section. 

Let us start with the bifurcation approximation of Eq.(2.18). Near the critical point 
B{x) -C A, the integrand can be safely replaced by a linearized expression, xB{x)/{x + 
B[xY) — > B[x) for dominant region of the integral x>M'^ with M = B{M^). Small 
but still important nonlinear effect from the region x < can be properly taken into 
account by introducing the infrared cutoff M in Eq.(2.18) except for some subtlety to 
be discussed later in the pure NJL limit. 

Under this approximation, the differential equation Eq.(2.19a) and UVBC Eq.(2.20b) 
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remain the same, while the Eq.(2.19b) and IRBC Eq. (2.20a) are changed into 



^ {xV{x)) = B{x) , 

T/(M2) = 0, S(M^) = M , 



(5.1a) 
(5.1b) 



respectively. Prom equation Eq.(5.1a) and Eq.(2.19a) we obtain the differential equation 

for V{x): 

/\ h Z^X \ 

^■^^ » V{x) = . (5.2) 



Before starting the full analysis of Eq.(2.18) with bifurcation approximation, we 
briefly describe the result in the l/Nf.^ hmit, in which the analytical calculation is 
greatly simplified without changing the qualitative feature. In this limit, the differential 
equation Eq.(5.2) and its boundary conditions are simplified: 



^ {xV{x)) + ^V{x) = 0, 
V{M^) = 0, B{M^) = M, 



B{x) 



d_ 
dx 



Ml 



[xVix)), 



-A. 



(5.3) 



The solution of Eq.(5.3) is given by 



V{x) 
B{x) 



M 
2/1 
M 
2/1 



-1/2+// 



/ X \-V2-M 



X \ 



(5.4a) 
(5.4b) 



where 2/i := \/l — 4A and M/A (= ^) is calculated 



as 



'1\ ^ ^ /M\^^ ^ 4/i-(l/2 + /x)% 
.C/J U; 4/i- (1/2 -/.)%' 



(5.5) 



with := M?/A^. Prom Eq.(5.5) we obtain the critical surface 



ih = (1/2 + /x)%. 



(5.6) 
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which is independent of Z^. Eq.(5.6) coincides with the critical hne of the gauged NJL 
model [27]. 

We next discuss the RG flow at the leading order. The function /i(A, M, F^, Aj? — 
A) of Eq.(4.1) can be read from Eq.(5.5): 

WA M F A) - ^^ (V2 + /x)^-(l/2-/i)^(M/Ar/- 

/i(A, M, F., A) - - 1 _ (M/A)4M ■ ^^-^^ 

From Eq.(3.9a)-(3.9b) the function Z<^(A, M, F^, A) is given by 

Z,(A, M. F„ A) = ^'■; ,(A.M.F!?AAr 
where the IR cutoff was introduced: 



^ Att'^Jm^ {x + B^{x)) 

since the mass function B{x) is determined only for x > in the bifurcation approxi- 
mation. Among these "RG" flows for various F^ and M, the "renormalized trajectory" 
possesses a distinguished property: The functions h and remain finite in the limit 
A/M oo (fixed point): 

h{A,M,F^,X) = h{A/M,l,F^/M,X) ^ h{oo,l,F^/M,X), (5.10a) 
Z^(A,M,F„A) = Z^(A/M,1,F,/M,A) ^Z^(oo,l,F,/M,A). (5.10b) 

Noting that ^(A^) for A/M oo, we can easily see from Eq.(5.7) and Eq.(5.8) 
that this condition can be satisfied only for 

F^ = F]{A' = oo), (5.11) 

otherwise — > ±oo for A/M — > oo. 
The mass function B{x) is independent of {h,Z^) which depends on A through 
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Eq.(5.7) and Eq.(5.8). Thus, it is also independent of the cutoff A and we can easily 
calculate the function F]{k^) := F|(oo) - F]{K^) by 



' Att^ Ja^ {x + B^{x))^ Att^ Ja^ X 

^ 7V,M^ (l/2 + /.H5/4-/./2) /MX^-^^ 

where we have made an approximation B{x > A^) = (M/2/i)(l/2 + //)(x/M^)~^/^+^. 

Plugging Eq.(5.11) and Eq.(5.12) into Eq.(5.7) and Eq.(5.8), we obtain the explicit 
form of the RG flow on the renormalized trajectory for M <^ A in the ~^ limit: 

MA,M,F,(oo),A) = ^(1/2 + /.^ {^ + JTfkr ' ^^-^^^^ 

Z,{A, M, Ffioo), X)h{A, M, F,(oo), A) = | (^4 - /^^/W + /^)^ ^ ^5_^3^^ 

where we estimated y(A^) ~ (M/2//)(M/A)^"^''. It should be emphasized that the 
disappearance of the "renormalized trajectory" at A = can be already observed at 
this stage: — > oo as A ^ 0. 

Now we return to the analysis of Eq.(2.18) for finite Nc- Rewriting Eq.(5.2) in terms 
of a new variable rj := Z^x/M"^ for v{rj): 



V{x) = r]-^'^+^'{l + r]f'^^''v{r]), (5.14) 

we obtain a differential equation for v: 

77(1 + ri)v"{ri) + (1 + 2/x + 2(1 + /X + v)ri)v' + ^(1 + 2/i)(l + 2v)v = 0, (5.15) 

where the prime denotes the differentiation with respect to 77, v' :— -^v, and fj, and u 
are defined by 



/i:= ^Vl-4A, u:= ^^Jl-4h/Z^Nc. (5.16) 
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The general solution of Eq.(5.15) is obtained as a linear combination of two hyper- 
geometric functions. The ratio of the two coefficients of the linear combination of the 
hypergeometric functions is determined from the IRBC Eq.(5.1b). Then we obtain the 
solution of Eq.(5.2): 

V{7i) = k[V47i)-V4v)], (5.17) 
where k is a normalization constant and V±{r)) is defined as 

ri \-^/2±M / N i/2±. F(i±//±// + ^'.i±//±//-^'.l±2//:-r7) 



(5.18) 



with TjM := Z^M^/M^ and 



uj:= -^l-4X-4h/Z^N,. (5.19) 



Then the mass function B{x) is obtained from Eq.(5.1a) as follows: 



(5.20) 



where V± are defined by 



V,i,) := (1/2 ± ,)/(i±/^^- + ^4±/^^--^.2±2;.;-,) 

The coefficient k is determined by the normalization condition B{M'^) — M. 

It is worth remarking that the asymptotic form of the solution for V{x) and B{x) 
are given as follows: 



V{ri) = Ciri--2+'' + C2V~-'-^, 



(5.22a) 
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B{x) = Ci(^ + /x)r/-HM + C'2(^-/i)r^'^ (5.22b) 



in the small region, and 



1/(77) = C[r]-^+'' + a.r]-^-'' , (5.23a) 
B{x) = C[{- + u;)7]-^+^ + C'^{- - u;)f]-^-'' (5.23b) 

in the large region, where Ci, C2, and C2 constants determined by the IRBC 
Eq.(5.1b) and the normalization condition S(M^) = M. 

From the UVBC Eq. (2.20b) we obtain the scahng relation, equi-^/ (:= A/M) surface, 
of the fermion dynamical mass for three parameters {h, A, Z^): 

\Aj [ Z^ + r^ ) \^ 

^ (4V>/ - (1/2 + n?) 17:^(1/2 + + V.iZ,./r,) z, 
(4Vr,-(l/2-/.)2)-(^(l/2-^)2 + T/_(V^,) V,^' 

where we have defined := M|/A^ and 

P(y) := ""If+^ + ' + ^i+^ + '-^'^ + '^'-^l . (5.25) 

When we take the limit 0, we obtain the scaling relation for the gauged NJL 

model[27]: 

(M\'^ _ {A + l/N;)h-r,{l/2 + ^,f 

U; (4 + l/Ar,)/i-r^(l/2-/x)2 ' ^ - > 

since P(0) = 1 and V±{Qi) — 0. This implies the critical line: 
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which yields the critical point of the pure NJL model in the limit A — > 0: 



4 + l/7Vc' 



The critical exponent um for the fermion mass defined by 



(5.28) 



(5.29) 



takes the continuously changing value (for A < 1/4): 



1 



i^M(A) = — 



4/x 2V1 -4A' 



(5.30) 



For further results on other critical exponents, see, e.g., Ref.[23]. 

In the vanishing gauge coupling limit, A — > 0, we obtain the scaling relations for the 
pure Yukawa model: 



7 M- I ^ 



(4r^ + AZ^ + llN,)h - (1/2 + u)Z^ 
+ 4Z^ + 1/N,)h- {1/2 -i^)Z^ 

2 tanh"^ 



l/2i/ 



exp 



ly /i 4h 



{4Z^+4r^+l/Nc)h-^Z^ 



4h 
Z^Nc 



exp 



exp 



AZ^ + Ar^-l/N, 



2mr + 2 tan"^ 



1 7 / 4fc 



(4Z^+4r^+l/Arc)/i-|20 



4/t 



- 1 



> 4h/N,) 



[Z^ = Ah/N,) 



(0 < < 4/l/iVe) 



(5.31) 

where n is an integer [n — 0,1,2,...). This scaling relation is not the same as the 
previous result [6] on the pure Yukawa model which does not include the tadpole con- 
tribution. Without tadpole we have the solution only for Z^ < Ah/N^. 
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In the — > oo limit, Eq.(5.31) reduces to 



M2 

A2" 



4/i- 1 
Ah 



— exp —2 tanh 




l + 8{h- 1/4) 



1 



(5.32) 



which imphes the critical point he — 1/4. On the other hand, the 



limit 



reproduces Eq.(5.28) for n — 0. Hence the critical line h — hc{Z^) extends from 
{h, A, Z^) = (r^/(4 + 1/Nc), 0, 0) to (1/4, 0, oo) in the A = plane. 

It is worth remarking that the scaling relations Eq.(5.26) and Eq.(5.31) do not cor- 
rectly reproduces the scaling relation of the pure NJL model by taking both limits 
— > and A — > [29], although the gauged Yukawa model reduces to the pure NJL 
model in this limit. This is because the bifurcation technique cannot be applied to the 
pure NJL model in this limit where the mass function has no dependence. Note 
that the above scaling relations obtained in the bifurcation technique give a meaningful 
result only in the neighborhood of the critical point M/A = which is not sensitive 
to the nonlinearity of the SD equation. Therefore the ^ limit of the scaling 
Eq.(5.26) and the A — > hmit of Eq.(5.31) cannot give the correct scahng relation, 
apart from the critical point and the critical exponent.^ 

Now let us consider the := A/M oo limit of the scaling relation Eq.(5.24), 
namely, the critical surface h — hc{X, Z^ in three-dimensional parameter space (/i. A, Z^) 
Since the left-hand side of Eq.(5.24) vanishes in the hmit for < A < Ac = 1/4 (// > 0), 
the critical surface obeys the equation 



^In the pure Yukawa model the critical exponent Vm defined by ^- = (Ji — hJ^Z^Y'^ changes 
continuously depending on the value in the case 0, while vm takes, irrespective of Z^, the 

mean-field value: vm = ^, as long as 7^ in any cutoff A, as was pointed out in Ref.[25]. In the 
latter case, all the critical exponents take their mean field values. Actually the chiral condensate takes 
the critical exponent 1/2, with which the numerical calculation shown in Fig. 2 is consistent. 



[4/i/r^ - (1/2 + ^if 



rcf, + Z^ 



(1/2 + + V+{Z^/r^)]P{Z^/r^) = 0, (5.33) 
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which imphes 

A ^ Z^ F{1 + + u + + n + u- 00,2 + 2fi; -Z^/r^) 
^2 ^' F{1 + n + u + + fi + u-uj,l + 2fi;-Z^/r^)_ 
^ F{^+ H + + + H + iy -uj,l + 2fi; -Z^/r^) ^ ^ 
F{1- iji-u + u;,^- ijL-u-u;,l-2iJ,;-Z^/r^) 

This is the exact form of the critical surface of the gauged Yukawa model. By solving 
this equation, we obtain the following critical surface: 

h = h.{X, Z,) = r,(l/2 + /x)^ . ^ .^^^^y^t V^c z, for (5.35) 

The above equation shows the deviation from the critical line Eq.(5.27) of the gauged 
NJL model {Z^ — 0) due to the presence of Z^. 

Next we consider the decay constant F-^. Using Eq.(2.13) and Eq.(2.22) with the IR 
cutoff M, we may rewrite Eq.(3.9a)-(3.9c) as 



with the function Ff{x) being defined by 

which is finite (for 0<A<Ac=l/4) from the asymptotic form Eqs.(5.22b),(5.23b) of 
the solution B{x). Fj{M'^) in the bifurcation plays the role of Fj°'^{:= F|(A^ = oo)) 
in the full nonlinear numerical analysis in Section 4. 

Let us now consider the curve (RG flow) on which both the decay constant and 
the fermion mass M := B{M'^) are fixed, in the three-dimensional parameter space 
{h, A, Z<^). There are three types of flows corresponding to the sign of F^ — F^{M'^). 
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For the case (II) > Fj{M'^) we have — > oo, because the condensation function 
V{h?) goes to zero as := A/M oo. In the case (I) < FjlJVP) we can have 
= 0. This is consistent with the numerical result. 
However, for F^ — Fj{M'^) we obtain a finite value for Z^, since the first term of 
Eq.(5.36) is identically zero and the second term converges according to the asymptotic 
form of B{x) and V{x), Eqs.(5.22a)-(5.23b). Thus the surface F^ = Ff{M^) yields 
renormalized trajectories. Then we calculate the fixed line, i.e., the intersection of 
the surface of the renormahzed trajectories F^ — Fj{M'^) and the critical surface 
Eq.(5.34). On the critical surface (A/M oo), the renormalized trajectories (Eq.(5.36) 
for F2 = F]{M^)) take the form 



'2N,hj-V{x)^ 
rl B{x) [Bjx) + \ + A) V{x) 



16/i V{x)[B{x) -V{x)] 



a;=A2— >c!o' 



(5.38) 



where the differentiation of B{x) and V{x) are rewritten by using the differential equa- 
tion Eqs.(5.1a),(2.19b). 

Substituting the UVBC (2.20b) into Eq.(5.38) we obtain a relation among the pa- 
rameters (/i. A, Z^): 

^\ Z^ + r^ J\ Z^ + r^ 2{Z^ + r^y 2 / 

+16AzJlV!±±l±VA^, + ,_l\ 0. (5.39) 

\ Z^ + r^ ) 

Note that Eq.(5.39) is not the renormalized trajectories unless the parameters (/i. A, Z^ 
lie on the critical surface, since A/M ^ oo is already taken in Eq.(5.38). By requiring 
both Eq.(5.39) and Eq.(5.34) to be satisfied simultaneously, we finally obtain the fixed 
line which is shown in Fig. 9. Specifically, we find that the fixed line extends from 
(/i. A, Z^) = (— ^g_)_4/jy , 1/4, 0) to (1/4,0, oo) when — 1. This is one of the main 
results of this paper. 
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Now we state the result on the fixed fine. The fixed fine on the critical surface 
separates = Ff{M^) surface into two regions, one of which extends to —>■ oo 
as A/M ^ oo and another to ^ for A/M > 1. The surface = Ff{M'^) 
may be identified with the surface of renormalized trajectories. So we can imagine the 
renormalized trajectories lying on the surface — Fj{M'^) terminate at the fixed line 
on the critical surface. 

This picture is consistent with the numerical result in Section 4. Thus the analytical 
study confirm the existence of the renormalized trajectories and the fixed line. It is 
again emphasize that — > cx) as A — > 0, i.e., no (nontrivial) fixed point in the pure 
Yukawa model. The fixed line is only revealed by presence of the gauge coupling. 

6 Conclusion and discussion 

In this paper we have investigated phase structure of the {SU{Nc)-, U{1)—) gauged 
Yukawa model with a global symmetry SU{2)l x SU{2)r, using the ladder SD equa- 
tion with "standing" gauge coupling, Eq.(2.7a), and the newly derived "generalized 
Pagels-Stokar (PS) formula" Eqs.(3.9a)-(3.9c). In the three-dimensional space of bare 
parameters (/i. A, Z^), we have obtained the ^/-constant surface (^/ := A/M, fermion 
correlation length), or the scahng relation (see Fig. 5, Eq.(5.24)). Then, as a hmit of 
^ cx) or the fermion dynamical mass M — > (A-fixcd), we obtained the critical 
surface (Fig. 5, Eq.(5.34)) which separates the spontaneous chiral symmetry breaking 
phase (M 7^ 0) and the chiral symmetric phase (M = 0). 

We have further obtained the "renormalization-group (RG) fiows" by requiring the 
NG boson decay constant F-^ as well as M be fixed as we change the cutoff A (Figs. 6, 7). 
In the three-dimensional space each RG fiow is a cross section of the equi-F^r surface 
shced by the A-constant plane. Among such RG fiows we have discovered the "renor- 
malized trajectories" which terminate on the critical surface (Figs. 6, 7). The intersec- 
tion of a set of renormalized trajectories with the critical surface may be identified 
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with a set of UV fixed points, the UV fixed fine (Fig.9), which exists only for nonzero 
gauge couphng (A > 0). In the Nc ^ oo hmit we obtained an exphcit analytical 
form of the renormalized trajectory as well as the fixed line, Eqs.(5.13a)-(5.13b). The 
analytical study has shown that the qualitative feature of the phase structure does 
not depend on Nc- The fixed line extends connecting (h,X,Z^) = (— Yg^i^yj^-, 1/4, 0) 
and (1/4, 0,oo). We thus have established existence of the fixed line in the gauged 
Yukawa model (A > 0) and non-existence of the (nontrivial) fixed point in the pure 
Yukawa model (A = 0) through the analytical study and the numerical one within the 
framework of the ladder SD equation. This strongly suggests that although the pure 
Yukawa model might be "trivial", the gauged Yukawa model may be a "nontrivial" 
(interacting) continuum field theory thanks to the presence of the gauge coupling. 

The above fixed fine should be compared with that in the gauged NJL model. The 
— >■ limit of the gauged Yukawa model is equivalent to the gauged NJL model with 
appropriate rescaling of the Yukawa coupling (see Eq.(5.27)). Now in the gauged NJL 
model with the standing gauge coupling, the fixed line is identified with the critical 
line itself [27, 29, 21], while the fixed line in the gauged Yukawa model deviates from 
that of the gauged NJL model (Z^ = 0) and drifts into the direction of > on the 
critical surface in the three-dimensional bare parameter space (Fig.9). Actually, the 
projection of the fixed line onto the — plane exactly coincides with the critical line 
of the gauged NJL model. This result is due to the propagating ("hopping") degree 
of freedom of the scalar field, since the corresponding scalar field in the gauged NJL 
model is merely an auxiliary one and has no kinetic term. 

It was pointed out [22] that, within the gauged NJL model, the line of M-constant 
does not exactly coincide with that of F^r-constant particularly in the weak gauge 
coupling region. Actually, M is a slowly decreasing function of the cutoff A along 
the F^-constant line as was demonstrated in the top quark condensate model[3, 2].^ 

^This is not a problem of the top quark condensate model as a phenomenological model with an 
explicit cutoff. This in fact was the very predictabihty of the top quark mass in that model. Actually, 
we can even take the formal limit A ^ oo of the model and obtain a finite ( "nontrivial" ) continuum 
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This suggests that we need (at least one) other relevant operators to obtain a sensible 
continuum field theory. Having introduced a degree of freedom, we can explicitly 
see how this problem is resolved in this enlarged coupling space. Let us discuss the 
small Z^{> 0) in the region (I): < Fj°, which is bounded between the surface of 
renormalized trajectories ending at the fixed line (F^r — Ff° surface) and the Z^f, — 
plane (see Fig. 7). If we specify a point on the RG flow and perform one step of RG 
(by decreasing the cutoff to the smaller ^f), the point moves to another point on the 
same flow with a larger Z^. Thus the requirement of both and M to be constant 
can only be met by the change in Z^p direction and the continuum theory is obtained 
at the fixed line with > 0. 

Here we should remark that our results (the existence of the critical surface, the fixed 
line and the renormalized trajectory) are based on the "standing" ansatz of the gauge 
coupling. The RG evolution of the scalar quartic coupling A,^ is also disregarded in 
the present paper. In principle, it is possible to take these effects into account in our 
analysis by fixing yet another "physical quantities" . The mass of the physical scalar 
boson m^^^ and the (infrared) scale of SU{Nc) gauge interaction Agauge would suit 
such a purpose. We then need to estimate the possible feedback from the running of 
those parameters to the restricted analysis here within the essentially two parameter 
space {h, Z^). 

As for the running of A<^, wc need to modify two aspects so as to take account of the 
feedback in the SD equation utilized in this paper. One is the inter-relation among 
{'ipip) and (a) which should be changed with respect to the quartic coupling (see 
Eq.(2.11)). This effect is, however, suppressed by (i/jiIj) / and thus we hope can be 
disregarded near the critical point as far as the phase transition is the second order. 
Another is the use of the A^ dependent dressed propagators of the scalar/pseudoscalar 
field in the rainbow graph of the scalar/pseudoscalar. This is also suppressed by 1/Nc 
and is expected to give only a tiny effect for large Nc- 



theory [22, 2], although such a hmit is pecuhar in the sense we just mentioned. 
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It is a rather delicate problem, however, to estimate the feedback from the (asymp- 
totically free) running of the gauge coupling A. It is straightforward to include such 
a running effect (at one-loop) into our analysis as we did in the gauged NJL model 
{Z^ — 0)[22] where almost all the RG flows tend to the pure NJL point (/i. A) = {he, 0) 
as A ^ oo; namely, the fixed line shrinks into a single point {h, A) = {he, 0), the UV 
fixed point, due to the one-loop running effect of the gauge coupling[22].^° Thus the 
fixed line of the gauged Yukawa model in {h, A, Z^) space may also shrink into A ^ for 
the running gauge coupling, since its projection onto the {h, A) plane is again expected 
to be that of the gauged NJL model as in the case of the standing gauge coupling 
discussed above. Now, this limit is rather subtle, since the fixed line extends into the 
region — > oo as A ^ 0, see Fig. 9. Thus even the qualitative structure of the RG 
seems to depend critically on the speed of running of the gauge coupling, or on the 
number of quark flavors for the QCD coupling. [22] Similarly, in the perturbative RG 
analysis of the standard model[30], the qualitative structure of the RG flow of the 
Yukawa coupling in fact depends on the speed of the running of the QCD coupling. 
This problem will be dealt with in the forthcoming paper. 

The phase structure obtained in this paper may be compared with results of the 
lattice theory. If we restrict the model to the pure Yukawa case (A = 0), our results 
on the phase structure are consistent with those of the lattice Yukawa model [31]. We 
here predict the phase structure of the gauged Yukawa model on the lattice using the 
parameterization which is similar to the lattice theory. After the transformation of the 
parameters {d: dimensionality of the space-time, a: lattice spacing): 

2dZ^ + a2M|' N,2dZ^ + o?Ml' ^ ' 

our Yukawa model (excluding the gauge part) can be cast into the lattice version with 



^°In Ref.[22] we found another UV fixed point (/i, A) = (0,0), tlie "pure QCD point", which has 
no correspondence to the fixed point of the model with the standing gauge coupling discussed in this 
paper. 
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the lattice action: 



S 




(6.2) 



X 



with i — 1, .., A^c and a — 0, ..,3 being indices of the SU{Nc) fundamental represen- 
tation and the vector representation of 0(4) = SU{2)l x SU{2)r, respectively. The 
phase structure of the gauged Yukawa model as well as the pure Yukawa model in this 
parameterization is shown in Fig. 10. Our claim on the existence of the fixed line in the 
presence of the gauge coupling is based on the nonperturbative analysis through the 
SD equation. Actually, such a remarkable feature does not appear by perturbatively 
including the effect of the gauge coupling [32]. In view of this there exists, to our 
knowledge, no available data from the lattice theory which are comparable with our 
result. Our finding in this paper will be sufficient to urge the lattice people to try 
to perform a full study of the gauged Yukawa model on the lattice, which will really 
confirm whether our claim is right or not. 

In this paper we have assumed > 0, since otherwise the vacuum would be unstable 
already in the scalar/pseudoscalar sector (there appear tachyon ghosts). In fact the 
ladder SD equation has a tachyon pole and cannot be Wick-rotated for < —r^. 
Furthermore the RG flow becomes singular at Z^ — — r^[l — {V^^c + 1 — l)/(4-/Vc)] 
(> —r^), which can be seen from the UVBC Eq. (2.20b) and the generalized Pagels- 
Stokar formula Eqs.(3.9a)-(3.9c). However, we can formally extend our RG flow to 
the negative value > Z^> — r^[l — {V^^c + 1 — l)/(4-/Vc)] in our analysis. Although 
this Z^ region might also be pathological, it is amusing to compare our RG flow with 
the lattice Yukawa model. By the above correspondence we can extend the RG flow 
of Fig. 10a (no gauge coupling) into the negative kh region, which indicates no sign 
of flxed point. This is consistent with the recent lattice analysis of the pure Yukawa 
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model. 

Finally, besides the RG analysis, we consider a physical application of the above 
information on the phase structure. Here we keep the UV cutoff A fixed. This is 
actually the case for the Yukawa-driven top mode model [6, 8] and the heavy scalar 
technicolor model[17]. First we discuss the criterion whether the symmetry breaking 
is "dynamical" or not. This is not so clear in the system having the elementary scalar 
from the onset, since it always mixes with the composite one. However, as we can 
see from our phase diagram, we can discriminate two distinct region separated by the 
renormalized trajectories: region (1) (F^^ < F^) and region (11) (F^ > F^). In region 
(I) we can always arrange the fermionic contribution F| to dominate the bosonic one 
F^ in the the NG boson decay constant F^, particularly for a small or a large A 
(see Fig. 8). This is contrasted with the region (II) where is relatively large and 
the bosonic contribution persists even for a large A. We may call such a fermionic 
dominance of the in region (I) the "dynamical", in contrast with the fermionic 
non-dominance in region (II) which may be called the "non-dynamical". Although 
the nomenclature is somewhat obscure, our definition is based on the phase diagram 
and hence is without ambiguity. This concept may be useful for analyzing the unified 
models based on the strong fermion dynamics involving elementary scalar field. 

Next we discuss the dependence of the fermion dynamical mass M on Z^. For a 
given value of F^, we have shown that M decreases as Z^ increases along the A- fixed 
line, see Fig. 8. This implies that M in the gauged Yukawa model {Z^ > 0) is always 
smaller than that in the gauged NJL model {Z^ — 0) for the same value, which is 
actually derived more generally as a direct consequence of our generalized PS formula 
Eqs.(3.9a)-(3.9c) as was explained in Section 3. This might appear rather peculiar 
from the general view point of a "generahzed NJL model" (without gauge coupling) of 
Hasenfratz et al.[15] presuming arbitrary higher dimensional operators. Such operators 
yield a priori an equal potentiality either to raise or lower the fermion dynamical mass 
compared with the NJL model. On the contrary, we have specified a model, the gauged 
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Yukawa model, for the higher dimensional operators instead of introducing arbitrary 
such operators. Therefore we have obtained a definite answer that the fermion mass 
in this model must be lower than in the gauged NJL model. If this conclusion in the 
SU{2)l X 5'C/(2)ij-invariant gauged Yukawa model persists in the SU{2)l x C/(l)y- 
invariant gauged Yukawa model, then the Yukawa-driven top quark condensate based 
on the latter model will give a smaller top quark mass prediction than the original 
one[3, 5] based on the SU{2)l x [/(l)y-invariant gauged NJL model. The analysis 
along this direction including the running effects of QCD coupling is under way. 
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Figure Captions 

Figure 1: 

Graphical representation for the NG boson decay constant Eq.(3.4). 
Figure 2: 

Critical line of the pure Yukawa model (A = 0). The critical line is drawn by the 
bold line and the solid lines denote the renormalizat ion-group flows. Each (upper) 
flow corresponds to a different (larger) value of F^. The dotted lines are the equi-^/ 
(correlation length) lines. 

Figure 3: 

Yukawa coupling dependence of the chiral condensate. The solid line, broken line, 
dashed line and dotted hne correspond to — 0.01, 0.1, 1, 10, respectively. 

Figure 4: 

Critical line (bold line) and equi-^/ lines for = (gauged NJL model). 
Figure 5: 

Equi-^/ (correlation length) surface in the space (/i, A, Z^. The framed surfaces from 
top to bottom corresponds to = A/M = exp(i:/2) for t = 1,3,5,7,15. The final 
surface is very close to the critical surface (^/ = oo). 

Figure 6: 

Renormalization- group flows for fixed F^, which are shown as the intersection of the 
equi-F^ surface with the fixed-A (gauge coupling ) plane, (a) F^/M"^ — 0.142, (b) 
F^/M^ = 0.086. 
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Figure 7: 

Critical line (bold line) and renormalization-group flows (solid line) for the gauged 
Yukawa model (A > 0). They are depicted in the fixed-A plane. Each (upper) flow 
corresponds to a different (larger) value of F^. The dotted lines denote the equi-^/ 
lines, (a) A = 0.05, (b) A = 0.15. 

Figure 8: 

The ratio F"^ / Fj along the renormalization-group flow in the fixed-A plane. The line 
passing through the origin denotes the fixed-cutoff A line, (a) A = 0.02, (b) A = 0.08. 

Figure 9: 

The mesh denoted by sold line shows the critical surface obtained from the analytical 
solution of the SD equation Eq.(2.19a) and Eq.(5.1a) in the limit of A/M — > oo. The 
fixed line on the critical surface is plotted by bold line. The projection of the fixed line 
to the Z^-h plane (a) and the \-Z^ plane (b) are plotted by sohd line. The equi-A hues 
(a) and cqui-^ lines (b) are shown by the lines corresponding to different values of A 
and h (see the upper right corner), respectively. 

Figure 10: 

Critical line (bold line) and renormalization-group flows (solid line) in the fixed-A plane. 
The parameters are translated by using Eq.(6.1), and we choose the lattice spacing as 
a^M| = 50 and Nc = 1. The dotted line denotes the equi-^/ line, (a) A = (pure 
Yukawa model), (b) A = 0.1. 
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